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ABSTRACT 


A three-dimensional, triply nested tropical cyclone forecast model 
was initialized from hand-analyzed synoptic-scale wind data. The diag- 
nostic phase forced the mass fields from the wind fields by use of a 
Suitable balance equation. Latent heat parameterization and frictional 
dissipation were omitted from the model to study movement primarily due 
to advective processes. For comparison, a three-level uniform coarse 
mesh grid model was initialized with the same real data. Time averaging 
of the pressure gradient terms of the momentum equations was incorpora- 
ted into each model in an attempt to increase the maximum time increment. 

Movement forecasts of Typhoon Irma to 48-hours with both dynamic 
models were compared with operational forecasts and the post-analysis 
"best" track. Very satisfactory movement prognoses were obtained for the 
three periods forecasted with the dynamic models. The models correctly 
forecast the passage of a 500-mb trough to the north of the typhoon, 
apparently preventing recurvature. For movement forecasts, it appears 
that the nested model can be initialized with linear interpolation of 
data from the coarse grid to the finer mesh grids. The time increment 
was increased through use of the pressure averaging in the uniform CMG 
model; however, ro time step increase was realized with the nested 


model due to lattice separation of the solution. 
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I. INTRODUCTION 


Accurate and timely prediction of tropical cyclone movement and in- 
tensity continues to be one of the most difficult realizations in the 
meteorological community. One reason for the absence of good numerical 
guidance is that the relatively large grid spacing used in today's 
numerical forecast models cannot resolve the tropical cyclone circulation. 
The application of a 381-kilometer (km) grid, such as that used by the 
Fleet Numerical Weather Central (FNWC), Monterey, California, not only 
enhances the phase error in short-wave systems, but also. lacks the reso- 
lution necessary to describe tropical cyclone energetics. 

Consequently, successful numerical simulation of mesoscale weather 
systems will require a finer grid resolution than that currently used in 
extratropical models. This extra resolution could be achieved by decreas- 
ing the grid increment of the mid-latitude model to properly define the 
mesoscale processes important to tropical cyclone development. However, 
necessary computer storage and model run time would be prohibitive. For 
example, for each reduction of the grid spacing by a factor of two, the 
computer time necessary to generate a forecast at the same temporal point 
is increased eight-fold. One alternative to a uniform fine mesh model 
is the use of multiple (nested) grid models. These models incorporate 
a smaller grid interval only in the immediate vicinity of the disturbance 
and utilize a coarser mesh grid elsewhere, An example is the three-level 
primitive-equation, triply nested model of Harrison (1973). Figure 1 


shows the grid nesting used in Harrison's model. Typically, the coarse 
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mesh grid (CMG) resolution is double that of the medium mesh grid (MMG). 
Likewise, the reduction of the grid spacing from the MMG to the fine 
mesh grid (FMG) is one-half. Clearly, the advantage of Harrison's (1973) 
nested model over one that employs a uniform fine mesh grid over its 
entire domain, is the considerable savings of model run time and com- 
pater Storage. 

Harrison's (1973) experiments were based upon analytically derived 
initial states. He found that heating in the FMG only, resulted in 
cyclonic development in the CMG. This development demonstrated that his 
two-way interaction between grids, will allow for momentum and heat trans- 
fer with only minimal effect on the solution. This important achieve- 
ment sets the stage for more sophisticated and diversified testing of 
the Harrison model. 

It is the purpose of this thesis to continue experiments on 
Harrison's model, directed toward the eventual development of an opera- 
tional tropical cyclone forecast model for the western North Pacific. 
The primary goal of this thesis is to test Harrison's (1973) model 
uSing operational real data obtained from the Fleet Weather Central/ 
Joint Typhoon Warning Center (FWC/JTWC) Guam, Mariana Islands. The 
Storm selected for this study was Typhoon Irma during the period 25 to 
27 November 1974. Throughout this time span, irma's location with res- 
pect to synoptic-scale features was SE that recurvature was a distinct 
possibility. Movement forecasts using the nested model and a uniform 
CMG model (unnested) will be examined to see if numerical guidance would 
have assisted the forecasters. It was anticipated that one of the pri- 
mary tests would involve the suitability of the western Pacific data 


network. 
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A real data initialization for a nested tropical cyclone model 
entails numerous general problem areas, the foremost of which is the 
procurement of data. Mathur (1974) produced a 96-hour movement forecast 
of Hurricane Isbel (1964) utilizing a nested version of Krishnamurti's 
(1973) model. His data was hand analyzed on the fine mesh grid scale, 
with the coarse mesh grid values extracted from the coincident fine mesh 
grid points. Such an initialization is ideal. That is, a model which 
forecasts mesoscale weather phenomena may be expected to perform best 
if real data were readily available on that scale. However, such is not 
the case in the western Pacific. A successful tropical cyclone model 
for this region must rely primarily on the synoptic scale for its data 
base. Then, this data will have to be interpolated to the finer mesh 
grids. If the storm is readily accessible, aircraft reconnaissance could 
support the initialization through collection of mesoscale data in the 
immediate vicinity of the cyclone. Assuming this data can be input to 
the model in a timely fashion, it would provide an excellent supplement 
to the synoptic-scale information. However, due to decreasing recon- 
naissance resources, the likelihood that sufficient mesoscale data can 
be collected is in question. Thus, initialization of Harrison's model 
will be attempted using hand analyzed synoptic-scale data extracted onto 
a 220-km grid. This CMG data will then be linearly interpolated to the 
finer mesh grids. 

A secondary objective of this thesis is to institute study of the 
feasibility of increasing the time step in Harrison's model. It was 
stated previously that the nested model's increased resolution near the 


storm vortex will necessitate increased model run time. This time 
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increase must somehow be reduced if the model is to be used operationally. 
Shuman (1971) and Brown and Campana (1971) found that by applying a 
weighed time average to the pressure gradient term of the momentum equa- 
tions, the maximum time step could be doubled. Clearly, successful 
application of this scheme to the nested model would help alleviate the 
increased time burden. Thus, Harrison's model was re-aligned to accept 
the time averaging scheme. Since the proposed scheme was not successful 
for nested versions, only a summary of the results will be presented in 


the Appendix. 
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II. PRIMITIVE EQUATION PREDICTION MODEL 


A. THE MODEL EQUATIONS 


The primitive-equation model after Harrison and Elsberry (1972) and 


Harrison (1973) was the basis for these experiments. 


are: 
OT OT 
Un ë x EC 
3t IDE EV M Ax + 3x 39 ; 
OT OT 
U E M o E yy 
"= L(v) fu M Jy + 3x + ay i 
90 
L -L (0) , 
991909 - D, 
ot B 1000 * 
QU) 2 ð u d v 
SSC Fi + 3 W] , 
ð$ 8 C d _p_,R/C 
dp p dp ‘1000 i 
where 
ZNOW NS d „VS d 
L(S) = M = cy? Í n» SL | 35 (ws) 


The model equations 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


L(S) represents the flux divergence of the general scalar quantity S. 


Other meteorological symbols used above can be found in the "List of 


Symbols." 


The linear computational stability criterion for two-dimensional 


equations governing simple wave motion is 


CAE 
— < ,7 
m 07 
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(7) 


where 


Q 
HI 


the phase velocity of the fastest gravity wave 


At 


HI 


time increment 


Le horizontal Grid increment 


Equation (7) dictates a maximum time step of 480 seconds in the nested 
CMG and the uniform CMG model. The MMG and FMG maximum time steps are 
240 and 120 seconds. An increase in the uniform CMG time increment to 
800 seconds was achieved by time averaging the pressure gradient term 
of the momentum equations (see Appendix). 

The initial step was forward in time. All subsequent iterations used 
the leapfrog scheme. Parameterization of heating and moisture was inten- 
tionally neglected from the system to emphasize storm movement primarily 


by advective processes. 


B THE GRID 

The grid utilized for this experiment followed that of Harrison 
(1973). It is a mercator secant projection true at 22.5°N. One can see 
from Fig. 1 that the grid extends from 100E to 154E and from the Equator 
to 38N. Figure l also illustrates the relative position of the nested 
model's MMG and FMG at the beginning of the integration on 25 November 
1974. The grid mesh for the uniform CMG and the CMG of the nested model 
was 220 km. Resolution on the nested model's MMG and FMG was 110 km and 
55 km, respectively. Figure 2, shows the vertical distribution of 
dependent parameters for both models. Although the variables are stag- 
gered in the vertical, horizontal space staggering was not utilized in 


these experiments. 
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CMG grid spacing is 220 km, MMG spacing is 110 km 


The triply nested as it appeared on 00 GMT 25 November 
and FMG increment is 55 km (after Harrison, 1973). 


Figure l. 
1974. 
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Fig. 2. Vertical distribution of dependent variables and 
pressure levels for the three-dimensional model 
(after Harrison, 1973). 
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C. BOUNDARY CONDITIONS 
Boundary conditions on the northern and southern walls for the uni- 
form CMG and nested CMG were the no-flux type after Elsberry and Harrison 


(T971). These conditions are 


om = d CH 
2o 

“541 M. = (8b) 

Er (8c) 


Such representation actually places the wall between the outer two grid 
rows on the poleward and equatorial sides. Mass is conserved along 
these boundaries since the v component is always zero at the wall, that 
is, the wall is a streamline. Thus, the 18 grid rows of extracted data 
were expanded to 20 rows with the inclusion of the no-flux boundaries. 
The east-west boundaries were made cyclic after Krishnamurti (1969). 
This procedure, applicable to the 1000-mb height field and the 850-, 500-, 
and 200-mb u,v fields, adds a buffer zone in the east-west direction to 
absorb the impact of forced continuity. The original data of 23 longi- 
tudinal columns was expanded to 28 columns by adding one column to the 
west side and four columns to the east side. The filling of this buffer 


is accomplished letting 


S . = S. (9a) 
2009 273 


, e S . (9b) 
Ee 21,3 


and then linearly interpolating to fill values in columns 25, 26 and 27. 


Here S represents any of the original data fields mentioned earlier. 
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Bucovenbpiesieveloerty at the upper boundary 1s equal to zero and is 
calculated at levels 5, 3 and 1, Fig. 2, through downward integration 
of Eq. (5). Optimum results were obtained when, in Eq. (4), o3 was 


replaced by the equivalent expression -wa . Further solution improvement 


was attained when w in Eq. (4) was represented by an averaged value from 


E = SET” RE SE Ger RER 200 MIA (10) 
PEZEBTERUSION 

Horizontal diffusion is a computational device that was used to con- 
trol small scale instability, to retard lattice solution separation and 
to remove high frequency noise from the model solution. It was not in- 
tended to represent parameterization of a physical atmospheric process, 
but rather to be utilized as a higher order smoothing mechanism. The 


diffusion components in the u momentum equation are 


A 324 
ne = D Ku Se (lla) 
AR SE 


Here K testne horizontal aiffusion coefficient. A similar represen: 
tation of diffusion was applied to Eq. (2). In the uniform CMG model, 
ei = 1 x aso. provided the best results. This coefficient was 
carefully chosen, because a value too large will smooth the meteorologi- 


cal modes as well as the spurious 2Ax and 3Ax waves. In the nested model, 


Eu = 4 x 107m /sec, was used on the MMG. The FMG value was 
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52 : , 
K — Ix 10m /sec. The difference by a multiple of four arises due to 
the dependence of the diffusion on the grid increment (see Eq. (lla)). 
No horizontal diffusion was used in the nested model CMG. Likewise, no 


Kertiecisdırftusion or frictional effects at the surface have been incor- 


porated into either model. 


E: EIENER 

Prior to output, the 1000-mb geopotential and wind velocity fields 
at all levels were smoothed twice with Shapiro's (1971) "ideal" 9-point 
filter. The advantage of this particular filter is that it completely 
removes 2Ax waves and most of the 3Ax modes. The filter is described 


generally by: 


dÉ 
— nn + 
Q 1 (2 (0 ) 4a. 


: Dot Sehe, 4; OCL. : 
14-15 zer, je gua asch 1,) 
(12) 


EE SE : IC: ; 


It was not employed during the integration process to avoid smoothing 


the meteorological waves. 
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III. NESTED GRID INTERACTION AND GRID MOVEMENT 


A. TWO-WAY INTERACTION VS ONE-WAY INTERACTION 

The primary feature that distinguishes nested grid models is the 
nature of the interaction between grid scales. Information can be 
passed between the different grid scales by either one-way interaction 
or two-way interaction (Phillips and Shukla, 1973). With one-way in- 
teraction, the larger scale solution determines the smaller scale 
motion, without being influenced by the changes occurring within the 
finer mesh grid. The iterative procedure is simple, since the grids 
are integrated separately using the same coding. This type of boundary 
best fits situations where the large scale physics are passive with 
respect to those on the smaller scale. That is, one-way interaction 
simply provides increased resolution with no feedback from the FMG to 
BE CHG 

It is well documented from previous studies, that an integral part 
of the tropical cyclone development process is the energy transfer 
from the small scale convective towers to the synoptic scale (convec- 
tive instability of the second kind (CISK)). If a numerical model is 
to Simulate the tropical cyclone energetics, it too should provide a 
means for energy transfer among the various scales of motion. Obviously, 
Since the smallest mesh length in the Harrison model is 30 nmi, energy 
transfer from the cloud scale to the mesoscale cannot be parameterized; 
however, the two-way interaction inherent in the model does allow for 


transfer of energy between the mesoscale (FMG) and synoptic scale (CMG). 
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Two-way interaction, such as that used by Harrison (1973), requires 
simultaneous integration on all grids. Figure 3 illustrates this inte- 
gration procedure. One can see that first the CMG is integrated one 
time step, then the MMG, then the FMG two steps, the MMG again and 
finally two more FMG iterations. This formal procedure is necessary 
since each reduction in the grid increment requires a corresponding 
change in the time step (Eq. (7)) and all grids must be at the same tem- 
poral location prior to the next CMG iteration. After each integration 
of the coarser mesh grid, the boundaries on the finer mesh grid are 
specified through temporal and spatial interpolation of the correspond- 
ing coarse grid values. At the conclusion of every FMG integration, 
the interior coincident grid points in the MMG grid are revised via a 
nine-point average, The same procedure is used to update the CMG with 
the MMG solution. Thus, the coarser grid provides the boundary values 
to the finer grid, which, upon completion of the finer mesh grid inte- 


gration, feeds back new values to the interior of the coarse grid. 


B. MOVEMENT OF THE NESTED GRIDS 

Successful movement of the finer mesh grids within their coarser 
counterparts is a necessary condition that must be fulfilled before 
long-term tropical cyclone forecasts can be generated with a nested 
model. The best solution on the FMG should be obtained if the storm be- 
ing monitored remains near the center of the grid. It can be surmised 
that close proximity of the storm to the grid boundary will lead to a 
deteriorating solution and possibly even computational instability. 
The longitudinal distance across the 17 x 17 point FMG with grid reso- 


INtionmemeso nmi is 510nmi. "A storm initially at the center of the grid 
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ae 


Figure 3. A schematic to illustrate the marching process in a nested 
grid model (Harrison, 1973). 
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Figure 4. An illustration of MMG and FMG movement. 
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will reach the west edge in one day, assuming a 10-knot average westward 
speed. Therefore, grid movement is an essential component of the 
Harrison model. 

When the minimum 1000-mb geopotential value corresponding to the 
approximate storm vortex moves to the seventh grid point from the western 
boundary of the FMG, the FMG grid is moved one MMG grid point to the west. 
From Fig. 4, one can see that the procedure entails the placement of the 
coincident MMG values onto their proper FMG positions. Then linear 
spatial interpolation is used to fill the remaining FMG voids. The two 
FMG eastern-most columns are disregarded and the storm center becomes 
relocated at the ninth FMG point from the west boundary. 

If, in turn, the FMG moves to within three grid ST the MMG west 
boundary, the MMG is moved one CMG point to the west. The new western 
MMG values are taken from the CMG and linearly interpolated in space to 
fill the remaining MMG vacancies. Referring to Fig. 4, movement of the 
MMG in the CMG also requires a redefining of the FMG location within the 
MMG. Identical provisions are provided for relocation of the FMG and 


MMG to the north, 
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IV. INITIALIZATION 


A. DATA INPUT 

Data for the initialization were supplied by the U. S. Fleet Weather 
Central/Joint Typhoon Warning Center, Guam, Mariana Islands (FWC/JTWC). 
The surface pressure and 850-, 500- and 200-mb isotachs and streamlines 
were hand-analyzed at 00 GMT 25, 26 and 27 November 1974. These analyses, 
plus the careful extraction of the data to the grid, provided a data base 
superior to objective analyses which might have been available. 

For convenience, the real atmospheric data at 500- and 200-mb were 
extracted to the nearest model predictive level (see Fig. 2). Subsequent 
reference to these levels will be according to the real data level, not 


the model's predictive level. 


B. TROPICAL VS EXTRATROPICAL 

In general, initialization of tropical models differs from their ex- 
tratropical counterparts. In a mid-latitude model, the wind is normally 
derived from a solution of the balance equation with an objectively ana- 
lyzed geopotential field (Haltiner, 1971). Such a procedure is not suit- 
able for a tropical model since the pressure gradient in the southern 
latitudes is very weak and the pressure field is not well defined. More- 
over, it can be shown that the mass field primarily tends to adjust to 
the wind field in the tropics during the early stages of the integration 
after initialization, especially in mesoscale models as found by Anthes 
(1974). Consequently, in the tropics, primary consideration must be given 


to the wind field when performing the initial diagnostic phase. 
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C. UNIFORM CMG DIAGNOSTIC PHASE 
The relative vorticity fields were obtained from the observed u and 


v components with 


eee UE 
C = MES (2) Ę 


r um M = EM )] (13) 


Subsequently, the streamfunction, U, was found by sequential over-relaxa- 


tion of the expression 


V EMG (14) 


uj =- 5, vj = ME (15) 


Solution of Eq. (14) requires specification of boundary values on the 
northern and southern peripheries. To remain consistent with the no-flux 
boundaries, described by Eqs. (8b) and (8c), constant streamfunction 
values were defined for the northernmost and southernmost two rows. 
Following Elsberry and Harrison (1971), zero was chosen for the southern 
ss and the northern values were calculated from 


Z u 


m 
Ve, cy Ay us) 


y 


north ` 


Here m is the mean zonal wind component averaged over the entire grid 


and Ay is the distance between the north and south walls. 


2% 





An appropriate balance equation can be derived through partial differ- 
entiation of Eq. (1) with respect to x and Eq. (2) with respect to y. 


Addition of these two equations, plus rearrangement of terms leads to 


u Lu" EV) 


V 
ee SM [= — nn Hm? = O a 


dt 30M y x un ces 
2 39 3 3^ 

-M up. 5 )-M ps + — (17) 
Y Ox e 


To begin the prognostic stage smoothly, the left hand side of Eq.(17) is 
set equal to zero. Surface diffusive effects were neglected.  Manipula- 


tion of Eq. (17) leads to a Poisson equation, 


2 R a L (u,,) a 2) Vu a ov ) f 
Vó=- I ]- FSWET — r D- nou 3, W (18) 


that is readily solved with sequential over-relaxation or other direct 
methods. In this case, the mean geopotential value on the south wall was 
derived hydrostatically (Eq. (6)) based upon a mean climatological tem- 
perature sounding. An average geopotential value for the north wall was 
found through integration of the v momentum equation (noting that I = 0) 


over the region, 





" = _ u, Xu.) 
$, 7 6, - (£ + C) + —— ay (19) 


Again, Ay is the north-south grid distance. To assure maintenance of 
the no-flux boundary conditions, Eq. (8a) was now applied at all levels. 
Having found the balanced geopotential fields from the upper level wind 


fields, only the 1000-mb geopotential field and the temperature fields 
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Bee yes 24, + and 6, Fig. 2, need to be specified to begin the 
prognostic phase. Since the initial wind fields are non-divergent, 
the vertical velocity is everywhere zero. 

The 1000-mb geopotential fielä was obtained hydrostatically from 
the 850-mb values. This requires that the 1000-mb and 850-mb potential 
temperature fields are defined. From Eq. (6), recall that the thick- 
ness between two levels in the atmosphere is proportional to the mean 
potential temperature of the layer. Constant potential temperatures for 
the 1000-mb north and south walls were obtained from climatology. A 
linear gradient at 1000 mb was forced to fit the boundary values. It 
was assumed that the 850-mb potential temperatures were everywhere 4°K 
warmer than the corresponding 1000-mb values. This assumption resulted 
in a linear meridional gradient in the 1000/850-mb layer mean potential 
temperature field. Finally, the 1000-mb geopotential field and the 
500-mb and 200-mb potential temperature distributions were obtained 


using Eq. (6). 


D. NESTED GRID MODEL DIAGNOSTICS 

The diagnostic procedure for the nested model CMG followed that of 
the uniform CMG model (see Fig. 5). A discussion of the diagnostic 
phase of the interior grids is now presented. Particular emphasis will 
be placed on the MMG and FMG initialization. 

l. Telescoping 

The divergent wind components on the MMG were interpolated 

linearly from the CMG values. The same interpolation scheme was applied 
to define the FMG wind distribution. Vorticity fields were calculated 


on both finer mesh grids. The streamfunction field was saved from the 
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previously relaxed CMG, and through linear interpolation, provided 
boundary values and an initial guess field for the solution of Eq. (14) 
on the MMG. Having calculated the non-divergent wind on the MMG, the 
FMG streamfunction initial guess field was formed, the solution obtained 
and the winds computed. 
Subsequently, the non-linear balanced equation was solved for 
the geopotential fields. Again the coarser grid solution provided the 
boundary values and initial guess fields for the finer mesh grid. The 
MMG and FMG 1000-mb and 850-mb potential temperature domains were inter- 
polated successively from the coarser grid. The 1000-mb geopotential 
field and the upper level potential temperature fields for the nested 
grids were derived hydrostatically. Finally, like the CHE, the MMG and 
FMG vertical velocity fields were equal to zero. 
2. malĘWeSSER Le lescoBing 

Elsberry and Harrison (1971) found that, in a limited region 
model, improper specification of the geopotential boundaries may lead to 
oscillations in the mass field. This "sloshing" results from a mismatch 
of the boundaries and the wind distribution. To eliminate the possible 
imbalance caused by the specification of the CMG mass and wind fields 
aS boundary values for the MMG and FMG, a procedure called "reverse" 
telescoping was introduced. 

Reverse telescoping makes use of an enlarged MMG and FMG during 
the diagnostic phase. Remember that optimum prognostic results would 
be expected if the model were initialized only from mesoscale data. 
The application of an expanded grid allows for greater usage of meso- 
scale data, which hopefully, will lead to minimization of the adjust- 


ment during the initial stages of the integration. 
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The foundation for the reverse telescoping is that during re- 
laxation of a Poisson equation, the effect imposed by the specified 
boundary values quickly diminishes toward the grid interior. This 
means that, although the total north-south geopotential gradient may 
not support the mean wind field between (see Eq. (19)), the interior 
grid will be consistent with the mean wind field implied by the forcing 
function. Thus if the MMG and FMG are expanded to a larger size, as 
NB the initialization within the interior normal 17 x 17 grid 
boundaries should represent an improved balanced state. MMG and FMG 
grid expansions to 21 x 21 points were used during the initialization. 
A still larger grid expansion may be expected to yield higher quality 
results; however, its use will probably be limited by the large increase 
in computer storage requirements. 

The application of reverse telescoping does not necessitate 
revision of the telescoping initialization described earlier, except 
to apply this procedure to the enlarged MMG and FMG. What is required, 
however, is an extension of the initialization in the reverse direction. 
After determining the FMG streamfunctions, or geopotential fields, new 
solutions for the MMG and CMG streamfunction and geopotential fields 
are calculated with the finer mesh grid values as interior boundaries. 
Re-application of the relaxation technique is necessary to ensure that 
the boundary conditions of the finer mesh grid correspond to the values 
in the coarser mesh grid coincident with the finer grid boundary. This 
preserves the two-way interaction between grids. 

To avoid repetition, the reverse phase will now be described as 


schematically illustrated in Fig. 7, assuming that the FMG streamfunction 
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Figure 6. Grid enlargement including depiction of the area relaxed 
during reverse telescoping. 





field has been solved. Subsequently, the inner 17 x 17 FMG streamfunc- 
tion values and non-divergent winds were extracted and placed at appro- 
priate positions on the MMG, The expanded MMG streamfunction grid was 
calculated within area A, Fig. 6. This ties the interior streamfunction 
boundary of the FMG with the remainder of the MMG field, and necessi- 
tates recalculation of the MMG non-divergent wind field in area A. 
Finally, the MMG streamfunction and non-divergent winds were extracted to 
their proper CMG positions. The CMG streamfunction field was determined 
within the perimeter area and the corresponding wind field was recomputed. 
Note, that the procedure required that the vorticity fields on all grids 
were retained from the telescoping phase of the calculations. 

Having finished the reverse E ders for the streamfunction, 
the method was similarly extended to the balance equation. Here the 
impact of the "reverse" telescoping should be more profound, because the 
interpolation of the data to the finer mesh grid was linear, and the 
balance equation contains several non-linear terms. Note also, that the 
telescoping phase for the geopotential fields relies on the winds obtained 
from the reverse telescoping of the streamfunction. 

Upon completion of the reverse telescoping only the 17 x 17 MMG 
and FMG geopotential and wind fields were retained for the prognostic 
stage. The combination of the telescoping and reverse telescoping phases, 


plus the remaining diagnostic steps, is summarized in Fig. 7. 
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V. RESULTS 


A. INITIAL AND 24-HR PROGNOSTIC FIELDS 

Example Of the initial streamfunction fields of the triply nested 
model at 00 GMT 25 November 1974 are presented in Figs. 8-10. Because 
the reverse telescoping procedure adds increased resolution primarily 
near the cyclone vortex, the uniform CMG model initial streamfunction 
distributions are nearly identical to those of the nested model. The 
massive low-level circulation is well depicted in Fig. 8. Since the 
models are initialized from hand-analyzed wind fields, there is good 
vertical consistency in the streamfunction analyses, especially near 
Typhoon Irma. Other important synoptic features include a moderately 
strong 500-mb short-wave trough northwest of Irma, and the subtropical 
high cells east and west of Irma, (Fig. 9). The 200-mb chart, Fig. 10, 
indicates a weak southeasterly flow over Irma, with no cyclonic vortex 
present at this stage. 

Having noted the similarity between the uniform CMG model non-diver- 
gent wind fields with those of the nested model, the effects of the 
reverse telescoping can be easily identified with the aid of Figs. 11 and 
12. One can see the effect of the increased resolution attained near 
the vortex in the reverse telescoped model (Fig. 12), and the Similarity 
of the fields in regions distant from the storm center. In this case 
the diagnosis of the mass field from the finer Scale wind fields intro- 
duces a closed 30-m contour that does not appear on the coarse grid field. 


The reversed telescoped MMG initial 1000-mb height distribution that 
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cmresponds to the CMG initial field in Fig. 12, is presented in Fig. 13. 
Increased definition of the cyclone vortex is clearly depicted here. 

Note that this grid contains the O-m and 30-m contours, supporting the 
fact that the higher resolution of the MMG allows for analysis of finer 
scale features. 

Prognostic fields from the triply nested model verifying at 00 GMT 
26 November are shown in Figs. 14-17. In many regards, the large scale 
wind fields are nearly identical to those prepared with the uniform coarse 
mesh model. The most significant 24-hour change, is the northeastward 
building of the 500-mb ridge to the west of Irma, as shown in Fig. 16. 
This occurs in response to the northerly flow behind the vacating mid- 
latitude trough observed northwest of Irma earlier, Be on the initial 
fields, it was expected by operational meteorologists that the advance 
of the 500-mb trough would produce a northward deflection and cause Irma 
to recurve into the westerlies. However, recurvature did not occur in 
the model (or in nature, see next Section) as the mid-latitude trough 
passed north of Irma. 

Careful observation of Figs. 15-17 will reveal some vertical separa- 
tion in the storm structure. That is, the 24-hour 500-mb forecast 
position is slightly north of the forecast center at lower levels. This 
vertical decoupling is caused by the shearing of the mean flow with height, 
the coarse resolution of the model, and the lack of deep convective heat- 
ing in the model. Greatest vertical decoupling was observed on 25 Novem- 
ber when the equatorward penetration of the westerlies was the greatest. 
Very little vertical separation was noted for the 24-hour forecast fields 


from 26 November, possibly suggesting that consistent steering flows 
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Figure 13. Nested model MMG initial 1000-mb height field at 00 GMT 
25 November 1974 (reverse telescoping). 
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existed in the vertical throughout this forecast period. Vertical separa- 
tion appeared to be slightly worse in the uniform CMG model, than in the 
nested model where the MMG and FMG tended to link the vertical structure 
near the storm vortex. One should note that the purpose of the study was 
to examine dynamic movement forecasts primarily due to advective processes, 
since vertical diffusion and latent heat release parameterization were 
omitted from the model. Numerical simulation of rain bands, the eye wall 
or even the intensity of the storm were not of concern in this experiment. 
Lastly, the degree of interaction between Irma and the mid-latitude 
features should be mentioned. Mathur's (1974) movement forecast of Hurri- 
cane Isbel essentially excluded interaction of the higher latitude forcing 
by imposing constant dependent variables on the northern boundary. In 
the Harrison model, the grid movement and two-way interaction between 
grids allows the storm to respond to the mid-latitude forcing. Further- 
more, the particular case chosen for this study represents a late-season 
storm when the extratropical influences are the greatest. The fact that 
the model successfully simulated the interaction in this advective case 
is encouraging; however, inclusion of latent heat parameterization and 
moisture to the model will be necessary before more conclusive comments 


can be made. 


Be TRACKS 
Forecasting of tropical cyclone movement is a very difficult task, 
and therefore requires that complete data resources and forecast tools 
are at the meteorologist's disposal. The primary data sources at FWC/JTWC 
are computer-derived tropospheric wind data, hand-analyzed charts at 


three atmospheric levels, detailed reconnaissance plots on sectional 
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charts, satellite data and radar data. Forecasting aids include analog 
forecasts, statistical forecasts with regression equations, extrapola- 
tion, Steering by geostrophic winds at 500 mb and 700 mb and steering by 
global band upper air fields. These objective techniques are further 
supplemented by varıous climatological publications. Thus, the prepara- 
tion of each tropical cyclone forecast by FWC/JTWC personnel entails 
much careful preparation and guidance from the current objective and 
statistical methods. The 24-, 48- and 72-hour typhoon movement and in- 
tensity forecasts are generally prepared at normal synoptic times, 00 
GMT plus every 6 hours (Annual Typhoon Report, 1974). 

Dynamic movement forecasts to 48 hours with the uniform CMG model 
and the nested model will be compared with FWC/JTWC ee tracks 
and the FWC/JTWC "best" track, as determined by a detailed post-season 
analysis. Fig. 18 depicts dynamic and operational track forecasts from 
OO GMT 25 November. FWC/JTWC forecast a recurvature track (Fig. 18, 
Track A), which 1s not surprising considering the strong break in the 
SUB ropical ridge at this time (Fig. 9]. Tracks B and D, Fig. 18, are 
those forecast by the uniform CMG model and the nested model, respec- 
tively. While the FWC/JTWC 24-hour forecast position is about 100 nmi 
Barth Of the "best" track position (Track C, Fig. 18), the uniform CMG 
and the nested forecast errors are considerably less, with the nested 
model prognostic position within 30 nmi of the verifying location. 
One can readily see that the 48-hour position forecast discrepancy 
created by the recurvature track of FWC/JTWC, might have been reduced 
considerably had numerical dynamic guidance been available. The uniform 
CMG model and nested model 48-hour forecast errors vary from about 40 nmi 


to near 60 nmi, with the better prognosis obtained with the uniform CMG 


model. 
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A comparison of the dynamic model prognoses, the FWC/JTWC "best" 
track and the FWC/JTWC forecast track from 00 GMT 26 November ZAWIE 
Paz ocnted in Fig. 19. Here deviation of about 40 nmi from the "best" 
track is observed in the initial position of the uniform CMG model and 
the nested model. This illustrates the difficulty of forcing the exact 
initial position onto the finite difference lattice when the FMG fields 
are simply linearly interpolated from the CMG analysis. The 24-hour 
forecast by FWC/JTWC is still a pre-recurvature track, although the 
slower movement indicates that a more westerly track was expected by 
FWC/JTWC forecasters, but was evidently not forecast to preserve con- 
tinuity with the previous prognosis (based on author's previous experi- 
ence). The nested model 500-mb streamfunction field at 00 GMT 26 Novem- 
ber 1974, shown in Fig. 20, reflects a 5-7 degree northward progression 
of the subtropical ridge in the 24 hours between the 25th and the 26th 
of November. During the same time period, Irma's 500-mb vortex moved 
slowly west-northwestward. Thus, the easterly steering component felt 
by Irma on the eguatorial side of the subtropical ridge should have 
increased. In fact, the "best" track, Track C shown in Fig. 19, re- 
veals that Irma moved slightly south of west during the next 48 hours 
and accelerated from an average 6-kt movement on the 25th to 10 kts on 
the 26th. The uniform CMG model forecast track very nearly parallels 
the "best" track; but is considerably slower. A slightly faster move- 
ment is forecast by the nested model, and consequently its 24-hour 
forecast position is somewhat closer to the "best" track position than 
either the uniform CMG model forecast position or the official FWC/JTWC 


forecast position. Errors for the uniform CMG 24-hour forecast and the 
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nested model 24-hour prognoses were about 140 nmi and 90 nmi, respective- 
ly. The advantage of the nested model over the uniform CMG model is 

lost by 48 hours, as the nested model produced a west-northwest track 

as opposed to the uniform CMG model west-southwest track that 

paralleled closely the "best" track. The uniform CMG 48-hour error was 
about 90 nmi while that of the nested model approached 135 nmi. 

The initial 500-mb streamfunction at 00 GMT 27 November is presented 
in Fig. 21. It shows continued slow northward progression of the sub- 
tropical high cells east and west of Irma. In general, however, the 
steering pattern appears to be very similar to that of the previous day. 
Dynamic guidance and the FWC/JTWC "best" track, as illustrated in Fig. 
22, both exhibit a west-southwesterly storm trajectory.  FWC/JTWC fore- 
cast a west-northwesterly track that resulted in 24- and 48-hour posi- 
tion errors of 66 nmi and 129 nmi, respectively (Annual Typhoon Report, 
1974). Position errors for the uniform CMG model were about 60 nmi 
at 24 hours and 160 nmi at 48 hours. The uniform CMG model 24-hour fore- 
cast position is slightly better than the 24-hour prognostic position 
of the nested model. It is interesting that most of the dynamic models' 
forecast errors arose due to an erroneous speed forecast, while the 
FWC/JTWC speed forecast was appreciably better and the direction fore- 
cast introduced most of the positioning error. 

In summary, the average 24- and 48-hour position errors for the 
uniform CMG model were approximately 85 nmi and 100 nmi. The same 
figures for the nested model were about 65 nmi and 160 nmi.  FWC/JTWC 


24- and 48-hour average values during this period were 115 nmi and 
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305 nmi; and for the 1974 season, 114 nmi and 218 nmi. Thus, it appears 
that the dynamic numerical model position forecasts during this period 
would have assisted operational forecasters. One should note, however, 
that the particular case studied here was selected to determine if the 
dynamic model would predict recurvature for this situation. It cannot 

be expected that the dynamic model would show such an improvement in 

every Situation. Continued testing of the nested model under varied move- 
ment situations, such as recurvature, is required before definite conclu- 


sions can be reached. 


O NELTECESCOPING VS REVERSE TELESCOPING 

Recall from the earlier chapter on initialization that linear inter- 
polation of dependent variables from the CMG to the finer mesh grids was 
adequate to specify the initial data fields prior to solution of the 
streamfunction and balance equations. Reverse telescoping was introduced 
to guarantee the ultimate balance between the mass and wind fields by 
solution on an enlarged grid. Relaxation in the reverse direction, that 
is, from the FMG, to the CMG, was then necessary to assure consistency in 
the grid interface values used in the model's two-way interaction. 

There are two benefits inherent in the reverse telescoping procedure. 
First, the interface boundary values are now determined by the finer 
mesh grid wind and height fields. As a result, the interface wind speeds 
in the finer mesh grids are characteristic of the higher resolution field 
and are not constrained by the coarser grid solution. During early steps 
of the integration, the interface boundary conditions are determined by 
the time tendency of the coarser mesh grid fields, but these tendencies 
are added to a finer mesh grid boundary value that is consistent with 
the interior field, Thus, the reverse telescoping will promote the two- 


wav Grid interaction, Second, relaxation on the enlarged grid reduces 
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the effect of the boundary values imposed by the CMG solution on the 
FMG. This allows for the best FMG dynamic balance between the mass and 
winds that is possible without using a fine mesh initialization every- 
where. Thus, the adjustment between the mass and wind distributions, 
during the first few iterations,due to dynamic imbalance in the initial 
fields, will be minimized. By avoiding the introduction of computa- 
tional modes early in the integration, contamination of the large-scale 
features is reduced. 

The reverse telescoping procedure contributes increased resolution 
to the MMG through the feedback of the FMG solution. This can be seen 
by comparing Figs. 23 and 24, the reverse telescoped and telescoped MMG 
initial 1000 mb height fields at 00 GMT 27 November. Note the expansion 
of the O-m and 30-m contours in the reverse telescoped field. The 
straight line segments of the -30-m isoline in Fig. 24 are due to the 
computer plotting package. Greater differences in the fields are not 
observed because non-linearity imposed during solution of the coarser 
mesh grid provides a fairly good boundary estimate to the finer mesh 
grid. More profound differences between the reverse telescoped and 
telescoped initial states may be realized in the initialization of 
asymmetric tropical cyclones, where non-linear gradients are observed 
and can only be represented by the data extraction to the coarser resolu- 
Cion grids: 

The improvement in the nested model prognostic fields when the 
initial state included reverse telescoping is not easily detectable in 
the large-scale fields. Relatively small changes were apparent in the 
forecast fields, but these were masked by filtering that proceeded the 


computer plotting. Furthermore, most computational modes introduced 
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Figure 23. Nested model MMG initial 1000-mb height field at 00 GMT 
27 November 1974 (reverse telescoping). 
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Nested model MMG initial 1000-mb height field at 00 GMT 
27 November 1974 (telescoping). 
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during the initialization can be seen only in the early stages of the 
integration and will be dampened by 24 hours. Lastly, the model's 
horizontal diffusion inhibits short-wave growth and may quickly stabi- 
lize gravity modes generated during the first few time steps. 

Two consequences of the reverse telescoping were readily observed 
in the initial data fields. First, the extraction of the winds from 
the finer mesh grid to the coarser grid resulted in an initial divergent 
component adjacent to the interface and within the coarser mesh grid 
area specified by the finer mesh grid solution. The resultant vertical 
velocities were, however, quickly dispersed early in the integration. 
Second, the feedback of data from the finer grids to the coarser mesh 
grids allows for increased compatibility between Bependent variables 
on the three grids. Since the FMG values are extracted to the MMG, and 
the MMG values to the CMG, higher resolution is passed from the FMG to 
the CMG. This increased consistency among the grids should result in 
a smoother start and a reduction of gravity modes during the initial 


iterative steps. 


D. NESTED MODEL-RELATED FEATURES 

A serious model-related problem was the lattice separation in fields, 
Since the dependent variables are carried at all grid points. This 
lattice separation increased throughout the integration and the "checker- 
board" character of the solutions would have led to computational insta- 
bility. To combat this separation, horizontal diffusion was used in 
the momentum equations (Eqs. (1) and (2)). The diffusion coefficients 
were larger than those necessary when analytical fields were used in the 


nested model (Harrison, 1973). The presence of horizontal diffusion 
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resulted in a gradual expansion of the maximum wind radius as time pro- 
gressed. This occurred on all grids, and represents a type of scale 
interaction. By 48 hours, the FMG cyclonic flow was greatly weakened, 
with nearly linear wind shear present in the horizontal. Once consider- 
able weakening in the zonal and meridional wind components was 
experienced, the FMG geopotential solution deteriorated. This explains 
why care must be taken when choosing the diffusion coefficients to be 
applied on the MMG and FMG. Inclusion of heating would have compensated 
the spurious dampening. However, it was not the intent of this research 
to test the various methods of latent heat parameterization which have 
been used in previous simulations of tropical cyclone development (see 
Eisberry, 1975). 

Another interesting computational aspect concerned the effects of 
grid movement. It was mentioned previously that when the tropical cyclone 
minimum 1000-mb height reached a position seven grid points from the FMG 
western or northern boundary, the FMG was moved one MMG point in the 
appropriate direction. Likewise, Similar provision was available to 
move the MMG one CMG point north or west when the FMG moved near the 
northern or western walls of the MMG. Each grid movement entails inter- 
polation of new geopotential values, potential temperatures and winds 
from the coarser mesh grid to the leading edge of the finer mesh grid. 
This linear interpolation of the wind components generates imbalances 
between the mass and wind fields near the leading edge of the grid 
boundary. As the wind and mass fields adjust, gravity modes will be 
propagated toward the center of the grid. Therefore, it seemed advisable 
to minimize the number of grid movements by placing the grids so that the 
storm center was initially located toward the southeastern portion of 


the FMG, and the FMG was located near the center of the MMG. 
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VI. CONCLUSIONS 


The three-dimensional triply nested model of Harrison (1973) was 
initialized using real hand-analyzed synoptic-scale data. Mass fields 
were forced by the winds by using a balance equation suitable for the 
tropics. Data on the finer mesh grids of the nested model were ob- 
tained through linear interpolation of the CMG field. The initializa- 
tion of the nested model incorporated "reverse" telescoping, a proce- 
dure that provides for the optimum balance between the finer mesh wind 
and mass fields through use of a grid expansion during the diagnostic 
phase. Initialization in the reverse direction, from the FMG to CMG, 
was required to specify consistent dependent variables on the grid 
interface. This last phase is necessary to assure a smooth start of 
the two-way grid interaction in the model. 

Latent heat and moisture parameterizations were neglected to allow 
for movement forecasts based primarily on advection. The intent of 
the study was to determine if dynamic guidance would have been of 
assistance in forecasting tropical cyclone movement, not to simulate 
the energy-related features of the tropical cyclone. Vertical diffu- 
sion was omitted from the model, but horizontal diffusion was required 
to control lattice solution separation that resulted because the dependent 
variables were carried at each grid point. 

Successful movement forecasts with the nested model and a uniform 
coarse grid model were obtained at 24 and 48 hours from 00 GMT 25, 26 


and 27 November 1974. During this period, Typhoon Irma was forecast 
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by FWC/JTWC personnel to recurve into the mid-latitude westerlies.  How- 
ever, both numerical models produced a westerly track, similar to the 
one actually followed by Irma. In all three cases the dynamic model 
could Ge been helpful to FWC/JTWC forecasters. The average dynamic 
model 24- and 48-hour position forecast errors for the three days were 
substantially lower than those of the operational forecasts. 

Thus, if one is interested only in storm movement, it appears that 
the nested model may be initialized with linearly interpolated data im- 
posed on the MMG and FMG, provided the mass and wind fields are properly 
balanced. Rather good storm movement was also simulated with the 
uniform CMG model for this case. This is especially important since the 
uniform CMG model incorporated the Shuman (1971) and Brown and Campana 
(1971) time averaging of the presSure-gradient terms. It allowed for 
a time step about 70% greater than the normal time step for a centered- 
in space and time finite difference scheme. 

Close proximity of the uniform CMG model forecast positions with 
those of the triply nested model reflects the diffusive spindown of the 
MMG and FMG winds and the dependence of cyclone movement on the large- 
scale advective flow. 

It is anticipated that future work on the model will entail incor- 
poration of a staggered horizontal grid to eliminate lattice separation, 
and thus reduce the need for horizontal diffusion. Heating parameteri- 
zation should be tested with the real data. But, most importantly, more 
case studies involving various intensity tropical cyclones under a 


number of types of steering flow should be completed. 
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APPENDIX: An Experiment Designed to Increase the Time Step of the 
Harrison Model 

A. INTRODUCTION 

The importance of a fast, efficient numerical model is manifest in 
the continuous research in this area. Faster machines, improved coding 
and computational techniques will eventually allow more sophisticated 
finite difference schemes to be utilized. An interesting computational 
scheme was introduced by Shuman (1971) and Brown and Campana (1971) and 
has been recently programmed for the National Meteorological Center 
(NMC) model. 

Shuman (1971) and Brown and Campana (1971) found that replacement 
of the pressure gradient term of the momentum equations with a weighted 
time averaged version permits doubling of the maximum time step. 


The averaging can be described by: 


] 7 023) Gi SR 
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where T is an arbitrary time level and "a" is a variable coefficient 
between O and 0.25. A similar form is used in the "v" momentum equation. 
When a = 0.25, theoretical results indicate a doubling of the time step 
can be achieved. Successively smaller time step increases can be 
realized as the coefficient a is reduced. When a = O, there exists no 
time averaging; that is, one reverts to use of the normal pressure 
gradient term. Experiments by Brown and Campana (1971) indicate that 
computational instability may arise if a = 0.25 is used. This insta- 
bility can be avoided by specifying a value of "a" slightly less than 


0.25 or through application of a suitable dampening mechanism. 
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B. PROCEDURE 

Since ay must be defined before Eq. (A~1) can be used, the 
solution order of the closed system was re-arranged to the following: 

1) lst Law of Thermodynamics 

2) 1000 mb boundary condition 

3) Hydrostatic equation 

4) Momentum equations 

5) Continuity equation 
For the condition, a - O, the rearrangement should yield the same 
answer as calculations using the normal equation sequence. This served 
as a check for programming errors. 

The time averaging was programmed successively into a two-dimensional 
(north/south cross-section) uniform grid model, a three-dimensional 
uniform CMG model, a two-dimensional nested model with a five-to-one 
grid reduction and the three-dimensional triply nested model.  Analytic 
initial states were imposed on all models and heating was applied on 
the finest grid scale of the two- and three-dimensional nested models 
as in Harrison and Elsberry (1972), and Harrison (1973). Model dependent 
averages of the w , in the o term of Eq. (4), were included. 


Diffusion, with various coefficients, was applied in the two-and three- 


dimensional nested models to control lattice solution separation. 


Ca T RESULTS 

Application of the time-averaging scheme to the two-and three- 
dimensional non-nested models was quite successful. In the three- 
dimensional uniform CMG model with grid spacing of 320 km, the maximum 


time increment allowed by Eq. (7) is approximately 750 sec. This 
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assumes that the fastest gravity modes move at 300 m/sec. The time 
averaged model solution closely approximated the Harrison model when 
the time increments for the two models were 1440 sec and 720 sec, 
respectively. All of the uniform CMG results in the results section 
of the main text were obtained from the time averaged model with a 
time step of 800 sec. This represents a 70% increase over the maximum 
time increment dictated by Eq. (7), and agrees well with Schoenstadt 
and Williams (1975), who found that the time advantage of the time 
averaging scheme is reduced to 70%-80% once a mean flow is added. 

Somewhat less satisfactory results were achieved when the time 
averaging scheme was inserted into the two- and three-dimensional nested 
models. Computational instability quickly destroyed solutions in both 
models when a time step greater than the maximum allowed by Eq. (7) 
was used. 

The two-dimensional nested model is characterized by only one grid 
reduction, a factor of five from the CMG to the FMG. When time averag- 
ing led to instability, the programming was checked by excluding the 
FMG solution portion from the iterative procedure, The instability 
appeared to be caused by lattice solution separation which was reduced 
by inclusion of horizontal diffusion when solving the momentum equations 
on the FMG. However, such strong diffusion was required that the 
diffusive time limit offset all advantage gained through the time averag- 
ing. To avoid application of horizontal qiEFusien to dampi the Nigh fres 
quency waves generated in the time averaged model, an alternative 
dampening mechanism was sought. 

At first, it was hoped that some type of energy absorbing scheme 


applied near the boundary in the FMG would absorb and prevent reflection 
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and growth of the gravity modes. A simple filtering scheme to reduce 
two grid length noise of the near boundary FMG v component was applied. 


This filter which used the three-point formula 


V, ie )/4 (A-2) 


+ 2v, + v. 
itl 1 _ 


1-1] 
provided no improvement. The second scheme, spongy boundaries, which 
allowed the FMG points adjacent to the interface to feel a percentage 

of the linear forcing of the CMG at the FMG boundary, also contributed 
little improvement. 

Since an analtyic initial state was used in the three-dimensional 
nested grid experiments, the forcing of the solution was provided by 
heating on the FMG. The long term solution of Harrison's model was 
nearly identical to that of the time averaged model when the same time 
step was used in both. This seems to suggest that the FMG heating func- 
tion forced the solution to an extent that the time averaging effect 
was eliminated. Furthermore, no time step increase beyond that 
dictated by Eq. (7) was achieved with the time averaging scheme in the 
triply nested model. Since the initiator of the computational insta- 
bility appears to be the lattice solution separation, the scheme will 
again be tested in a staggered grid nested model. Additional work 
appears to be necessary before use of the time-averaging scheme is 


employed in a nested grid model. 
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